3.4.25 \(\int \frac {\text {sech}^3(c+d x)}{a+b \sinh ^2(c+d x)} \, dx\) [325]
Optimal. Leaf size=92 \[ \frac {(a-3 b) \text {ArcTan}(\sinh (c+d x))}{2 (a-b)^2 d}+\frac {b^{3/2} \text {ArcTan}\left (\frac {\sqrt {b} \sinh (c+d x)}{\sqrt {a}}\right )}{\sqrt {a} (a-b)^2 d}+\frac {\text {sech}(c+d x) \tanh (c+d x)}{2 (a-b) d} \]
[Out]
1/2*(a-3*b)*arctan(sinh(d*x+c))/(a-b)^2/d+b^(3/2)*arctan(sinh(d*x+c)*b^(1/2)/a^(1/2))/(a-b)^2/d/a^(1/2)+1/2*se
ch(d*x+c)*tanh(d*x+c)/(a-b)/d
________________________________________________________________________________________
Rubi [A]
time = 0.07, antiderivative size = 92, normalized size of antiderivative = 1.00, number of steps
used = 5, number of rules used = 5, integrand size = 23, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.217, Rules used = {3269, 425, 536,
209, 211} \begin {gather*} \frac {b^{3/2} \text {ArcTan}\left (\frac {\sqrt {b} \sinh (c+d x)}{\sqrt {a}}\right )}{\sqrt {a} d (a-b)^2}+\frac {(a-3 b) \text {ArcTan}(\sinh (c+d x))}{2 d (a-b)^2}+\frac {\tanh (c+d x) \text {sech}(c+d x)}{2 d (a-b)} \end {gather*}
Antiderivative was successfully verified.
[In]
Int[Sech[c + d*x]^3/(a + b*Sinh[c + d*x]^2),x]
[Out]
((a - 3*b)*ArcTan[Sinh[c + d*x]])/(2*(a - b)^2*d) + (b^(3/2)*ArcTan[(Sqrt[b]*Sinh[c + d*x])/Sqrt[a]])/(Sqrt[a]
*(a - b)^2*d) + (Sech[c + d*x]*Tanh[c + d*x])/(2*(a - b)*d)
Rule 209
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[b, 2]))*ArcTan[Rt[b, 2]*(x/Rt[a, 2])], x] /;
FreeQ[{a, b}, x] && PosQ[a/b] && (GtQ[a, 0] || GtQ[b, 0])
Rule 211
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]
Rule 425
Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] :> Simp[(-b)*x*(a + b*x^n)^(p + 1)*
((c + d*x^n)^(q + 1)/(a*n*(p + 1)*(b*c - a*d))), x] + Dist[1/(a*n*(p + 1)*(b*c - a*d)), Int[(a + b*x^n)^(p + 1
)*(c + d*x^n)^q*Simp[b*c + n*(p + 1)*(b*c - a*d) + d*b*(n*(p + q + 2) + 1)*x^n, x], x], x] /; FreeQ[{a, b, c,
d, n, q}, x] && NeQ[b*c - a*d, 0] && LtQ[p, -1] && !( !IntegerQ[p] && IntegerQ[q] && LtQ[q, -1]) && IntBinomi
alQ[a, b, c, d, n, p, q, x]
Rule 536
Int[((e_) + (f_.)*(x_)^(n_))/(((a_) + (b_.)*(x_)^(n_))*((c_) + (d_.)*(x_)^(n_))), x_Symbol] :> Dist[(b*e - a*f
)/(b*c - a*d), Int[1/(a + b*x^n), x], x] - Dist[(d*e - c*f)/(b*c - a*d), Int[1/(c + d*x^n), x], x] /; FreeQ[{a
, b, c, d, e, f, n}, x]
Rule 3269
Int[cos[(e_.) + (f_.)*(x_)]^(m_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]^2)^(p_.), x_Symbol] :> With[{ff = Free
Factors[Sin[e + f*x], x]}, Dist[ff/f, Subst[Int[(1 - ff^2*x^2)^((m - 1)/2)*(a + b*ff^2*x^2)^p, x], x, Sin[e +
f*x]/ff], x]] /; FreeQ[{a, b, e, f, p}, x] && IntegerQ[(m - 1)/2]
Rubi steps
\begin {align*} \int \frac {\text {sech}^3(c+d x)}{a+b \sinh ^2(c+d x)} \, dx &=\frac {\text {Subst}\left (\int \frac {1}{\left (1+x^2\right )^2 \left (a+b x^2\right )} \, dx,x,\sinh (c+d x)\right )}{d}\\ &=\frac {\text {sech}(c+d x) \tanh (c+d x)}{2 (a-b) d}-\frac {\text {Subst}\left (\int \frac {-a+2 b-b x^2}{\left (1+x^2\right ) \left (a+b x^2\right )} \, dx,x,\sinh (c+d x)\right )}{2 (a-b) d}\\ &=\frac {\text {sech}(c+d x) \tanh (c+d x)}{2 (a-b) d}+\frac {(a-3 b) \text {Subst}\left (\int \frac {1}{1+x^2} \, dx,x,\sinh (c+d x)\right )}{2 (a-b)^2 d}+\frac {b^2 \text {Subst}\left (\int \frac {1}{a+b x^2} \, dx,x,\sinh (c+d x)\right )}{(a-b)^2 d}\\ &=\frac {(a-3 b) \tan ^{-1}(\sinh (c+d x))}{2 (a-b)^2 d}+\frac {b^{3/2} \tan ^{-1}\left (\frac {\sqrt {b} \sinh (c+d x)}{\sqrt {a}}\right )}{\sqrt {a} (a-b)^2 d}+\frac {\text {sech}(c+d x) \tanh (c+d x)}{2 (a-b) d}\\ \end {align*}
________________________________________________________________________________________
Mathematica [A]
time = 0.17, size = 91, normalized size = 0.99 \begin {gather*} \frac {-2 b^{3/2} \text {ArcTan}\left (\frac {\sqrt {a} \text {csch}(c+d x)}{\sqrt {b}}\right )+2 \sqrt {a} (a-3 b) \text {ArcTan}\left (\tanh \left (\frac {1}{2} (c+d x)\right )\right )+\sqrt {a} (a-b) \text {sech}(c+d x) \tanh (c+d x)}{2 \sqrt {a} (a-b)^2 d} \end {gather*}
Antiderivative was successfully verified.
[In]
Integrate[Sech[c + d*x]^3/(a + b*Sinh[c + d*x]^2),x]
[Out]
(-2*b^(3/2)*ArcTan[(Sqrt[a]*Csch[c + d*x])/Sqrt[b]] + 2*Sqrt[a]*(a - 3*b)*ArcTan[Tanh[(c + d*x)/2]] + Sqrt[a]*
(a - b)*Sech[c + d*x]*Tanh[c + d*x])/(2*Sqrt[a]*(a - b)^2*d)
________________________________________________________________________________________
Maple [B] Leaf count of result is larger than twice the leaf count of optimal. \(271\) vs.
\(2(80)=160\).
time = 2.02, size = 272, normalized size = 2.96
| | |
method |
result |
size |
| | |
risch |
\(\frac {{\mathrm e}^{d x +c} \left ({\mathrm e}^{2 d x +2 c}-1\right )}{d \left (a -b \right ) \left (1+{\mathrm e}^{2 d x +2 c}\right )^{2}}+\frac {i \ln \left ({\mathrm e}^{d x +c}+i\right ) a}{2 \left (a -b \right )^{2} d}-\frac {3 i \ln \left ({\mathrm e}^{d x +c}+i\right ) b}{2 \left (a -b \right )^{2} d}-\frac {i \ln \left ({\mathrm e}^{d x +c}-i\right ) a}{2 \left (a -b \right )^{2} d}+\frac {3 i \ln \left ({\mathrm e}^{d x +c}-i\right ) b}{2 \left (a -b \right )^{2} d}+\frac {\sqrt {-a b}\, b \ln \left ({\mathrm e}^{2 d x +2 c}+\frac {2 \sqrt {-a b}\, {\mathrm e}^{d x +c}}{b}-1\right )}{2 a \left (a -b \right )^{2} d}-\frac {\sqrt {-a b}\, b \ln \left ({\mathrm e}^{2 d x +2 c}-\frac {2 \sqrt {-a b}\, {\mathrm e}^{d x +c}}{b}-1\right )}{2 a \left (a -b \right )^{2} d}\) |
\(241\) |
derivativedivides |
\(\frac {\frac {\frac {2 \left (\left (\frac {b}{2}-\frac {a}{2}\right ) \left (\tanh ^{3}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )+\left (-\frac {b}{2}+\frac {a}{2}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{\left (\tanh ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )+1\right )^{2}}+\left (a -3 b \right ) \arctan \left (\tanh \left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{\left (a -b \right )^{2}}+\frac {2 b^{2} a \left (\frac {\left (-a +\sqrt {-b \left (a -b \right )}+b \right ) \arctan \left (\frac {a \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )}{\sqrt {\left (2 \sqrt {-b \left (a -b \right )}-a +2 b \right ) a}}\right )}{2 a \sqrt {-b \left (a -b \right )}\, \sqrt {\left (2 \sqrt {-b \left (a -b \right )}-a +2 b \right ) a}}-\frac {\left (a +\sqrt {-b \left (a -b \right )}-b \right ) \arctanh \left (\frac {a \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )}{\sqrt {\left (2 \sqrt {-b \left (a -b \right )}+a -2 b \right ) a}}\right )}{2 a \sqrt {-b \left (a -b \right )}\, \sqrt {\left (2 \sqrt {-b \left (a -b \right )}+a -2 b \right ) a}}\right )}{\left (a -b \right )^{2}}}{d}\) |
\(272\) |
default |
\(\frac {\frac {\frac {2 \left (\left (\frac {b}{2}-\frac {a}{2}\right ) \left (\tanh ^{3}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )+\left (-\frac {b}{2}+\frac {a}{2}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{\left (\tanh ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )+1\right )^{2}}+\left (a -3 b \right ) \arctan \left (\tanh \left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{\left (a -b \right )^{2}}+\frac {2 b^{2} a \left (\frac {\left (-a +\sqrt {-b \left (a -b \right )}+b \right ) \arctan \left (\frac {a \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )}{\sqrt {\left (2 \sqrt {-b \left (a -b \right )}-a +2 b \right ) a}}\right )}{2 a \sqrt {-b \left (a -b \right )}\, \sqrt {\left (2 \sqrt {-b \left (a -b \right )}-a +2 b \right ) a}}-\frac {\left (a +\sqrt {-b \left (a -b \right )}-b \right ) \arctanh \left (\frac {a \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )}{\sqrt {\left (2 \sqrt {-b \left (a -b \right )}+a -2 b \right ) a}}\right )}{2 a \sqrt {-b \left (a -b \right )}\, \sqrt {\left (2 \sqrt {-b \left (a -b \right )}+a -2 b \right ) a}}\right )}{\left (a -b \right )^{2}}}{d}\) |
\(272\) |
| | |
|
|
|
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(sech(d*x+c)^3/(a+b*sinh(d*x+c)^2),x,method=_RETURNVERBOSE)
[Out]
1/d*(2/(a-b)^2*(((1/2*b-1/2*a)*tanh(1/2*d*x+1/2*c)^3+(-1/2*b+1/2*a)*tanh(1/2*d*x+1/2*c))/(tanh(1/2*d*x+1/2*c)^
2+1)^2+1/2*(a-3*b)*arctan(tanh(1/2*d*x+1/2*c)))+2*b^2/(a-b)^2*a*(1/2*(-a+(-b*(a-b))^(1/2)+b)/a/(-b*(a-b))^(1/2
)/((2*(-b*(a-b))^(1/2)-a+2*b)*a)^(1/2)*arctan(a*tanh(1/2*d*x+1/2*c)/((2*(-b*(a-b))^(1/2)-a+2*b)*a)^(1/2))-1/2*
(a+(-b*(a-b))^(1/2)-b)/a/(-b*(a-b))^(1/2)/((2*(-b*(a-b))^(1/2)+a-2*b)*a)^(1/2)*arctanh(a*tanh(1/2*d*x+1/2*c)/(
(2*(-b*(a-b))^(1/2)+a-2*b)*a)^(1/2))))
________________________________________________________________________________________
Maxima [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Failed to integrate} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(sech(d*x+c)^3/(a+b*sinh(d*x+c)^2),x, algorithm="maxima")
[Out]
(a*e^c - 3*b*e^c)*arctan(e^(d*x + c))*e^(-c)/(a^2*d - 2*a*b*d + b^2*d) + (e^(3*d*x + 3*c) - e^(d*x + c))/(a*d
- b*d + (a*d*e^(4*c) - b*d*e^(4*c))*e^(4*d*x) + 2*(a*d*e^(2*c) - b*d*e^(2*c))*e^(2*d*x)) + 8*integrate(1/4*(b^
2*e^(3*d*x + 3*c) + b^2*e^(d*x + c))/(a^2*b - 2*a*b^2 + b^3 + (a^2*b*e^(4*c) - 2*a*b^2*e^(4*c) + b^3*e^(4*c))*
e^(4*d*x) + 2*(2*a^3*e^(2*c) - 5*a^2*b*e^(2*c) + 4*a*b^2*e^(2*c) - b^3*e^(2*c))*e^(2*d*x)), x)
________________________________________________________________________________________
Fricas [B] Leaf count of result is larger than twice the leaf count of optimal. 772 vs.
\(2 (80) = 160\).
time = 0.48, size = 1644, normalized size = 17.87 \begin {gather*} \text {Too large to display} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(sech(d*x+c)^3/(a+b*sinh(d*x+c)^2),x, algorithm="fricas")
[Out]
[1/2*(2*(a - b)*cosh(d*x + c)^3 + 6*(a - b)*cosh(d*x + c)*sinh(d*x + c)^2 + 2*(a - b)*sinh(d*x + c)^3 + (b*cos
h(d*x + c)^4 + 4*b*cosh(d*x + c)*sinh(d*x + c)^3 + b*sinh(d*x + c)^4 + 2*b*cosh(d*x + c)^2 + 2*(3*b*cosh(d*x +
c)^2 + b)*sinh(d*x + c)^2 + 4*(b*cosh(d*x + c)^3 + b*cosh(d*x + c))*sinh(d*x + c) + b)*sqrt(-b/a)*log((b*cosh
(d*x + c)^4 + 4*b*cosh(d*x + c)*sinh(d*x + c)^3 + b*sinh(d*x + c)^4 - 2*(2*a + b)*cosh(d*x + c)^2 + 2*(3*b*cos
h(d*x + c)^2 - 2*a - b)*sinh(d*x + c)^2 + 4*(b*cosh(d*x + c)^3 - (2*a + b)*cosh(d*x + c))*sinh(d*x + c) + 4*(a
*cosh(d*x + c)^3 + 3*a*cosh(d*x + c)*sinh(d*x + c)^2 + a*sinh(d*x + c)^3 - a*cosh(d*x + c) + (3*a*cosh(d*x + c
)^2 - a)*sinh(d*x + c))*sqrt(-b/a) + b)/(b*cosh(d*x + c)^4 + 4*b*cosh(d*x + c)*sinh(d*x + c)^3 + b*sinh(d*x +
c)^4 + 2*(2*a - b)*cosh(d*x + c)^2 + 2*(3*b*cosh(d*x + c)^2 + 2*a - b)*sinh(d*x + c)^2 + 4*(b*cosh(d*x + c)^3
+ (2*a - b)*cosh(d*x + c))*sinh(d*x + c) + b)) + 2*((a - 3*b)*cosh(d*x + c)^4 + 4*(a - 3*b)*cosh(d*x + c)*sinh
(d*x + c)^3 + (a - 3*b)*sinh(d*x + c)^4 + 2*(a - 3*b)*cosh(d*x + c)^2 + 2*(3*(a - 3*b)*cosh(d*x + c)^2 + a - 3
*b)*sinh(d*x + c)^2 + 4*((a - 3*b)*cosh(d*x + c)^3 + (a - 3*b)*cosh(d*x + c))*sinh(d*x + c) + a - 3*b)*arctan(
cosh(d*x + c) + sinh(d*x + c)) - 2*(a - b)*cosh(d*x + c) + 2*(3*(a - b)*cosh(d*x + c)^2 - a + b)*sinh(d*x + c)
)/((a^2 - 2*a*b + b^2)*d*cosh(d*x + c)^4 + 4*(a^2 - 2*a*b + b^2)*d*cosh(d*x + c)*sinh(d*x + c)^3 + (a^2 - 2*a*
b + b^2)*d*sinh(d*x + c)^4 + 2*(a^2 - 2*a*b + b^2)*d*cosh(d*x + c)^2 + 2*(3*(a^2 - 2*a*b + b^2)*d*cosh(d*x + c
)^2 + (a^2 - 2*a*b + b^2)*d)*sinh(d*x + c)^2 + (a^2 - 2*a*b + b^2)*d + 4*((a^2 - 2*a*b + b^2)*d*cosh(d*x + c)^
3 + (a^2 - 2*a*b + b^2)*d*cosh(d*x + c))*sinh(d*x + c)), ((a - b)*cosh(d*x + c)^3 + 3*(a - b)*cosh(d*x + c)*si
nh(d*x + c)^2 + (a - b)*sinh(d*x + c)^3 + (b*cosh(d*x + c)^4 + 4*b*cosh(d*x + c)*sinh(d*x + c)^3 + b*sinh(d*x
+ c)^4 + 2*b*cosh(d*x + c)^2 + 2*(3*b*cosh(d*x + c)^2 + b)*sinh(d*x + c)^2 + 4*(b*cosh(d*x + c)^3 + b*cosh(d*x
+ c))*sinh(d*x + c) + b)*sqrt(b/a)*arctan(1/2*sqrt(b/a)*(cosh(d*x + c) + sinh(d*x + c))) + (b*cosh(d*x + c)^4
+ 4*b*cosh(d*x + c)*sinh(d*x + c)^3 + b*sinh(d*x + c)^4 + 2*b*cosh(d*x + c)^2 + 2*(3*b*cosh(d*x + c)^2 + b)*s
inh(d*x + c)^2 + 4*(b*cosh(d*x + c)^3 + b*cosh(d*x + c))*sinh(d*x + c) + b)*sqrt(b/a)*arctan(1/2*(b*cosh(d*x +
c)^3 + 3*b*cosh(d*x + c)*sinh(d*x + c)^2 + b*sinh(d*x + c)^3 + (4*a - b)*cosh(d*x + c) + (3*b*cosh(d*x + c)^2
+ 4*a - b)*sinh(d*x + c))*sqrt(b/a)/b) + ((a - 3*b)*cosh(d*x + c)^4 + 4*(a - 3*b)*cosh(d*x + c)*sinh(d*x + c)
^3 + (a - 3*b)*sinh(d*x + c)^4 + 2*(a - 3*b)*cosh(d*x + c)^2 + 2*(3*(a - 3*b)*cosh(d*x + c)^2 + a - 3*b)*sinh(
d*x + c)^2 + 4*((a - 3*b)*cosh(d*x + c)^3 + (a - 3*b)*cosh(d*x + c))*sinh(d*x + c) + a - 3*b)*arctan(cosh(d*x
+ c) + sinh(d*x + c)) - (a - b)*cosh(d*x + c) + (3*(a - b)*cosh(d*x + c)^2 - a + b)*sinh(d*x + c))/((a^2 - 2*a
*b + b^2)*d*cosh(d*x + c)^4 + 4*(a^2 - 2*a*b + b^2)*d*cosh(d*x + c)*sinh(d*x + c)^3 + (a^2 - 2*a*b + b^2)*d*si
nh(d*x + c)^4 + 2*(a^2 - 2*a*b + b^2)*d*cosh(d*x + c)^2 + 2*(3*(a^2 - 2*a*b + b^2)*d*cosh(d*x + c)^2 + (a^2 -
2*a*b + b^2)*d)*sinh(d*x + c)^2 + (a^2 - 2*a*b + b^2)*d + 4*((a^2 - 2*a*b + b^2)*d*cosh(d*x + c)^3 + (a^2 - 2*
a*b + b^2)*d*cosh(d*x + c))*sinh(d*x + c))]
________________________________________________________________________________________
Sympy [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \int \frac {\operatorname {sech}^{3}{\left (c + d x \right )}}{a + b \sinh ^{2}{\left (c + d x \right )}}\, dx \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(sech(d*x+c)**3/(a+b*sinh(d*x+c)**2),x)
[Out]
Integral(sech(c + d*x)**3/(a + b*sinh(c + d*x)**2), x)
________________________________________________________________________________________
Giac [F(-2)]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Exception raised: TypeError} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(sech(d*x+c)^3/(a+b*sinh(d*x+c)^2),x, algorithm="giac")
[Out]
Exception raised: TypeError >> An error occurred running a Giac command:INPUT:sage2:=int(sage0,sageVARx):;OUTP
UT:Warning, need to choose a branch for the root of a polynomial with parameters. This might be wrong.The choi
ce was done
________________________________________________________________________________________
Mupad [B]
time = 6.14, size = 2797, normalized size = 30.40 \begin {gather*} \text {Too large to display} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(1/(cosh(c + d*x)^3*(a + b*sinh(c + d*x)^2)),x)
[Out]
exp(c + d*x)/((exp(2*c + 2*d*x) + 1)*(a*d - b*d)) - (2*exp(c + d*x))/((a*d - b*d)*(2*exp(2*c + 2*d*x) + exp(4*
c + 4*d*x) + 1)) + ((2*atan((b^2*exp(d*x)*exp(c)*(a*d^2*(a - b)^4)^(1/2))/(2*a*d*(a - b)^2*(b^3)^(1/2))) - 2*a
tan((exp(d*x)*exp(c)*((64*(20*a^3*d*(b^3)^(5/2) - 232*a^6*d*(b^3)^(3/2) + 2*a^9*d*(b^3)^(1/2) + 2*a*b^5*d*(b^3
)^(3/2) + 10*a^5*b*d*(b^3)^(3/2) - 22*a^8*b*d*(b^3)^(1/2) - 10*a^2*b^4*d*(b^3)^(3/2) - 20*a^4*b^2*d*(b^3)^(3/2
) - 18*a^2*b^7*d*(b^3)^(1/2) + 102*a^3*b^6*d*(b^3)^(1/2) - 242*a^4*b^5*d*(b^3)^(1/2) + 310*a^5*b^4*d*(b^3)^(1/
2) + 98*a^7*b^2*d*(b^3)^(1/2)))/(a*b^4*(a - b)^5*(a*b - a^2)*(a^2 - 2*a*b + b^2)*(a*d^2*(a - b)^4)^(1/2)*(3*a*
b^2 - 3*a^2*b + a^3 - b^3)*(9*a*b^2 - 6*a^2*b + a^3 - b^3)*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2
+ 6*a^3*b^2*d^2)^(1/2)) - (32*(b^8*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) +
36*a^2*b^6*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) - 47*a^3*b^5*(a^5*d^2 +
a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) + 30*a^4*b^4*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2
- 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) - 9*a^5*b^3*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3
*b^2*d^2)^(1/2) + a^6*b^2*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) - 12*a*b^7
*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2)))/(a^2*b^2*d*(a - b)^7*(a*b - a^2)*
(b^3)^(1/2)*(a^2 - 2*a*b + b^2)*(3*a*b^2 - 3*a^2*b + a^3 - b^3)*(9*a*b^2 - 6*a^2*b + a^3 - b^3)*(a^5*d^2 + a*b
^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2))) + (32*exp(3*c)*exp(3*d*x)*(b^8*(a^5*d^2 + a*b^4*
d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) + 36*a^2*b^6*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a
^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) - 47*a^3*b^5*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*
d^2)^(1/2) + 30*a^4*b^4*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) - 9*a^5*b^3*
(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) + a^6*b^2*(a^5*d^2 + a*b^4*d^2 - 4*a
^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2) - 12*a*b^7*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2
+ 6*a^3*b^2*d^2)^(1/2)))/(a^2*b^2*d*(a - b)^7*(a*b - a^2)*(b^3)^(1/2)*(a^2 - 2*a*b + b^2)*(3*a*b^2 - 3*a^2*b +
a^3 - b^3)*(9*a*b^2 - 6*a^2*b + a^3 - b^3)*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2
)^(1/2)))*((a^2*b^10*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2))/64 - (a^3*b^9*
(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2))/8 + (7*a^4*b^8*(a^5*d^2 + a*b^4*d^2
- 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2))/16 - (7*a^5*b^7*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*
a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2))/8 + (35*a^6*b^6*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3
*b^2*d^2)^(1/2))/32 - (7*a^7*b^5*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2))/8
+ (7*a^8*b^4*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2))/16 - (a^9*b^3*(a^5*d^2
+ a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2))/8 + (a^10*b^2*(a^5*d^2 + a*b^4*d^2 - 4*a^4*
b*d^2 - 4*a^2*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2))/64)))*(b^3)^(1/2))/(2*(a^5*d^2 + a*b^4*d^2 - 4*a^4*b*d^2 - 4*a^2
*b^3*d^2 + 6*a^3*b^2*d^2)^(1/2)) + (atan((exp(d*x)*exp(c)*(a^7*(a^4*d^2 + b^4*d^2 - 4*a*b^3*d^2 - 4*a^3*b*d^2
+ 6*a^2*b^2*d^2)^(1/2) - 3*b^7*(a^4*d^2 + b^4*d^2 - 4*a*b^3*d^2 - 4*a^3*b*d^2 + 6*a^2*b^2*d^2)^(1/2) + 55*a*b^
6*(a^4*d^2 + b^4*d^2 - 4*a*b^3*d^2 - 4*a^3*b*d^2 + 6*a^2*b^2*d^2)^(1/2) - 15*a^6*b*(a^4*d^2 + b^4*d^2 - 4*a*b^
3*d^2 - 4*a^3*b*d^2 + 6*a^2*b^2*d^2)^(1/2) - 297*a^2*b^5*(a^4*d^2 + b^4*d^2 - 4*a*b^3*d^2 - 4*a^3*b*d^2 + 6*a^
2*b^2*d^2)^(1/2) + 423*a^3*b^4*(a^4*d^2 + b^4*d^2 - 4*a*b^3*d^2 - 4*a^3*b*d^2 + 6*a^2*b^2*d^2)^(1/2) - 272*a^4
*b^3*(a^4*d^2 + b^4*d^2 - 4*a*b^3*d^2 - 4*a^3*b*d^2 + 6*a^2*b^2*d^2)^(1/2) + 90*a^5*b^2*(a^4*d^2 + b^4*d^2 - 4
*a*b^3*d^2 - 4*a^3*b*d^2 + 6*a^2*b^2*d^2)^(1/2)))/(a^8*d*(a^2 - 6*a*b + 9*b^2)^(1/2) + b^8*d*(a^2 - 6*a*b + 9*
b^2)^(1/2) + 130*a^2*b^6*d*(a^2 - 6*a*b + 9*b^2)^(1/2) - 314*a^3*b^5*d*(a^2 - 6*a*b + 9*b^2)^(1/2) + 367*a^4*b
^4*d*(a^2 - 6*a*b + 9*b^2)^(1/2) - 230*a^5*b^3*d*(a^2 - 6*a*b + 9*b^2)^(1/2) + 79*a^6*b^2*d*(a^2 - 6*a*b + 9*b
^2)^(1/2) - 20*a*b^7*d*(a^2 - 6*a*b + 9*b^2)^(1/2) - 14*a^7*b*d*(a^2 - 6*a*b + 9*b^2)^(1/2)))*(a^2 - 6*a*b + 9
*b^2)^(1/2))/(a^4*d^2 + b^4*d^2 - 4*a*b^3*d^2 - 4*a^3*b*d^2 + 6*a^2*b^2*d^2)^(1/2)
________________________________________________________________________________________